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Abstract 

There is growing evidence that some long gamma-ray bursts (GRBs) arise from the core 
collapse of massive stars, and thus it is inevitable that the environments of these GRBs are pre- 
burst stellar winds or dense media. We studied, for the first time, the wind model for afterglows 
based on the Blandford-McKee self-similar solution of a relativistic shock, and suggested that 
GRB 970616 is an interactor with a stellar wind. We also proposed a dense medium model for 
some afterglows, e.g., the steepening in the light curve of the R-band afterglow of GRB 990123 
may be caused by the adiabatic shock which has evolved from an ultrarelativistic phase to a 
nonrelativistic phase in a dense medium. We further discussed the dense medium model in more 
details, and investigated the effects of synchrotron self absorption and energy injection. A shock 
in a dense medium becomes nonrelativistic rapidly after a short relativistic phase. The afterglow 
from the shock at the nonrelativistic stage decays more rapidly than at the relativistic stage. 
Since some models for GRB energy sources predicted that a strongly magnetic millisecond pulsar 
may be born during GRB formation, we discussed the effect of such a pulsar on the evolution 
of the nonrelativistic shock through magnetic dipole radiation. We found that in the pulsar 
energy injection case, the dense medium model fits very well all the observational data of GRB 
980519. Recently, we combined the dense medium model with the pulsar energy injection effect 
to provide a good fit to the optical afterglow data of GRB 000301C. 

1 Introduction 

In the standard afterglow shock model (for a review see [27,38]), a gamma-ray burst (GRB) 
afterglow is usually believed to be produced by synchrotron radiation or inverse Compton scat- 
tering in an ultrarelativistic shock wave expanding in a homogeneous medium. As more and 
more ambient matter is swept up, the shock gradually decelerates while the emission from such 
a shock fades down, dominating at the beginning in X-rays and progressively at optical to radio 
energy band. The standard model is based on four basic assumptions: (1) the total energy of 
the shock is released impulsively before its formation; (2) the medium swept up by the shock is 
homogeneous and its density (n) is the one of the interstellar medium ~ 1 cm -3 ; (3) the electron 
and magnetic field energy fractions of the shocked medium and the index (p) in the accelerated 
electrons' power-law distribution are constant during the whole evolution stage; and (4) the 
shock is spherical. 

Each of these assumptions has been varied to discuss why some observed afterglows deviate 
from that expected by the standard afterglow model. For example, the R-band light curve of 
GRB 970508 afterglow peaks around two days after the burst, and there is a rather rapid rise 
before the peak which is followed by a long power-law decay. There are two models explaining 
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this special feature: (i) It was envisioned [28] that a postburst fireball may contain shells with 
a continuous distribution of Lorentz factors. As the external forward shock sweeps up ambient 
matter and decelerates, internal shells will catch up with the shock and supply energy into 
it. A detailed calculation shows that this model can explain well this special feature [26]. (ii) 
We considered continuous energy injection from a strongly magnetized millisecond pulsar into 
the shock through magnetic dipole radiation [8]. This model can also account for well the 
observations. It is very clear that these models don't use basic assumption (1). 

There are several models in the literature that discuss the effect of inhomogeneous media 
on afterglows [9,23,3,4], dropping the second assumption. Generally, anna r~ k (k > 0) medium 
is expected to steepen an afterglow's temporal decay. We studied, for the first time, the wind 
model for afterglows based on the Blandford-McKee self-similar solution [1] of a relativistic 
adiabatic shock, and suggested that GRB 970616 is an interactor with a stellar wind of n oc r~ 2 
[9]. It was found [3] that a Wolf-Rayet star wind likely leads to an n oc r~ 2 medium, and thus 
if GRB 980519 resulted from the explosion of such a massive star, subsequent evolution of a 
relativistic shock in this medium is consistent with the steep decay in the R-band light curve 
of the afterglow from this burst. Another way of dropping the second assumption is that the 
density of an ambient medium is invoked to be as high asn~ 10 6 cm -3 . The temporal decay 
of the R-band afterglow of GRB 990123 has been detected to steepen about 2.5 days after this 
burst [21,2,16]. We proposed a plausible model in which a shock expanding in a dense medium 
has evolved from a relativistic phase to a nonrelativistic phase [11]. We found that this model fits 
well the observational data if the medium density is about 3 x 10 6 cm~ 3 . We further suggested 
that such a medium could be a supernova or supranova or hypernova ejecta. 

In basic assumption (3), the electron and magnetic field energy fractions of the shocked 
medium may not be varied during whole evolution, as argued in [34], where all the observational 
data including both the prompt optical flash and the afterglow of GRB 990123 were analyzed. 

The steepening in the light curves of the afterglows of some bursts may also be due 
to lateral spreading of a jet, as analyzed in [29,31] when the jet expands in a homogeneous 
interstellar medium (ISM). This in fact drops basic assumption (4). However, numerical studies 
of [24,20,36,37] show that the break of the light curve is weaker and smoother than the one 
analytically predicted when the light travel effects related to the lateral size of the jet and 
a realistic expression of the lateral expansion speed are taken into account. In the case of a 
jet expanding in a wind, the calculated light curve is even much weaker and smoother than the 
ananlytical one [22,17]. We recently calculated light curves for GRB afterglows when anisotropic 
jets (dE/dQ oc 9~ k ) expand both in the interstellar medium and in the wind medium [5]. We 
found that in each type of medium, one break appears in the late-time afterglow light curve 
for small k but becomes weaker and smoother as k increases. When k > 2, the break seems to 
disappear but the afterglow decays rapidly. Thus, we expect that the emission from expanding, 
highly anisotropic jets provides a plausible explanation for some rapidly fading afteglows whose 
light curves seem to have no break. 

We discussed the dense medium model in more details [12], by taking into account both 
the synchrotron self-absorption effect in the shocked medium and the energy injection effect of 
[8,10]. Recently, we combined the dense medium model with the pulsar energy injection effect 
to provide a good fit to the optical afterglow data of GRB 000301C [13]. Here we want to give 
a brief review of some of our studies on GRB afterglows. 
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2 Shock Evolution 



It is well known that the evolution of a partially radiative shock depends on both the efficiency 
with which the shock transfers its bulk kinetic energy to electrons and magnetic fields and on 
the efficiency with which the electrons radiate their energy. In 1998, we proposed, for the first 
time, a unified model for dynamical evolution of a partially radiative shock [6]. This model 
is not only valid during the whole evolution stage including the Sedov phase for an adiabatic 
shock, but also can describe well an adiabatic shock as well as a highly radiatve shock. This 
model was later re-investigated and referred to as a generic one in [19]. For simplicity, we here 
assume that a relativistic shock expanding in a dense medium is adiabatic. This assumption is 
correct particularly for a low electron energy density fraction in the shocked medium [6,7]. The 
Blandford-McKee self-similar solution [1] gives the Lorentz factor of an adiabatic relativistic 
shock, 



1 

7= 4 



17£ (l + z) 



1/ 



lE 5 1 2 / 8 n 5 - 1 /V /8 [(l + ,)/2] 3 / 8 , (1) 



where Eq = E^ x 10 52 ergs is the total isotropic energy, n§ = n/10 5 cm -3 , i® is the observer's 
time since the gamma-ray trigger in units of 1 day, z is the the redshift of the source generating 
this shock, and m p is the proton mass. We assume 7 = 1 when % = tj,. This implies 

n 5 = E 52 t^[(l + z)/2} 3 . (2) 

For i® > tb, the shock will be in a nonrelativistic phase. In the following we will discuss the 
spectrum and light curve during the non-relativistic phase. 

As usual, only synchrotron radiation from the shock is considered. To analyze the spec- 
trum and light curve, one needs to know three crucial frequencies: the synchrotron peak fre- 
quency (u m ), the cooling frequency (i/ c ), and the self-absorption frequency {y a )- We assume a 
power law distribution of the electrons accelerated by the shock: dn' e /d"f e oc 7~ p for j e > 7 em , 
where 7 e is the electron Lorentz factor and 7 em = 610e e (7 — 1) is the minimum Lorentz factor. 
We further assume that e e and e b are the electron and magnetic energy density fractions of the 
shocked medium respectively. The v m is the characteristic synchrotron frequency of an electron 
with Lorentz factor of 7 em , while the v c is the characteristic synchrotron frequency of an electron 
which cools on the dynamical age of the shock. According to Sari et al. [32], we have derived the 
synchrotron peak frequency, the cooling frequency and the synchotron self-absorption frequency, 
measured in the observer's frame [12]. They are correct for the whole evolution stage. 

Now we give the spectrum and light curve of the afterglow during the non-relativistic 
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phase. First, for the case without energy injection, the shock velocity decays as oc i e and 
thus we have 



Fu = < 



W^)- (p - 1)/2 ^ m oc v -to-W t f-*M™ i£ Va <u<u e (3) 

{v c /v m )-^)/\u/v c )-^F Urn oc u-p/H^ 3p)/2 if v > u c . 



We easily see that for high-frequency radiation the temporal decay index a = (21 — 15p)/10 for 
emission from slow-cooling electrons or a = (4 — 3p)/2 for emission from fast-cooling electrons. 
If p w 2.8, then a ps —2.1 or —2.2. Comparing this with the relativistic result, we conclude that 
the afterglow decay steepens at the nonrelativistic stage. 
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Some models for GRB energy sources (for a brief review see [10]) predict that during the 
formation of an ultrarelativistic fireball required by GRB, a strongly magnetized millisecond 
pulsar will be born. If so, the pulsar will continuously input its rotational energy into the for- 
ward shock of the postburst fireball through magnetic dipole radiation because electromagnetic 
waves radiated by the pulsar will be absorbed in the shocked medium [8,10]. Since an initially 
ultrarelativistic shock discussed in [12] rapidly becomes nonrelativistic in a dense medium, we 
next investigate the evolution of a nonrelativistic adiabatic shock with energy injection from a 
pulsar. The total energy of the shock is the sum of the initial energy and the energy which the 
shock has obtained from the pulsar: 

E + / Ldt e = E tot oc v 2 r 3 , (4) 
Jo 

where L is the stellar spindown power oc (1 + i^/T) -2 (T is the initial spindown time scale). 
The term on the right-hand side is consistent with the Sedov solution. Please note that L can 
be thought of as a constant for t® < T, while L decays as oc t^ 2 for 3> T. Because of this 
feature, we easily integrate the second term on the left-hand side of equation (17). We now 
define a time at which the shock has obtained energy ~ Eq from the pulsar, t c = Eq/L, and 
assume t c <C T. The evolution of the afterglow from such a shock can be divided into three 
stages. 

Stage (i): i® <C t c , viz., the second term on the left-hand side of equation (4) can be 
neglected. The evolution of the afterglow is the same as in the above case without any energy 
injection. 

Stage (ii): for T > i® S> t c , the term Eq in equation (4) can be neglected. At this stage, 

—2/5 

the shock's velocity v oc i ffi , we have derived the spectrum and light curve of the afterglow 
[12] 

(z y a/z y m)~ ( ' P ~ 1 " ) / 2 (^ /v a ) b l 2 F Vrn OC V if v < v a 

{v/Vm)-^ 1)/2 F„ m OC ^-b-D^-W jf Va < v K Vc (5) 

(u c /^ m )-^ 2 (^c)- p/2 F Vm « v-vl 2 tl v if v > v c . 

It can be seen that for high-frequency radiation the temporal decay index a = (12— 5p)/5 ~ —0.4 
for emission from slow-cooling electrons or a = 2 — p ~ —0.8 for emission from fast-cooling 
electrons if p s=s 2.8. This shows that the afterglow decay may significantly flatten due to the 
effect of the pulsar. 

Stage (iii): for t® S> T, the power of the pulsar due to magnetic dipole radiation rapidly 
decreases as L oc , and the evolution of the shock is hardly affected by the stellar radiation. 
Thus, the evolution of the afterglow at this stage will be the same as in the above case without 
any energy injection. 

In summary, as an adiabatic shock expands in a dense medium from an ultrarelativistic 
phase to a nonrelativistic phase, the decay of radiation from such a shock will steepen, subse- 
quently may flatten if a strongly magnetic millisecond pulsar continuously inputs its rotational 
energy into the shock through magnetic dipole radiation, and finally the decay will steepen again 
due to disappearance of the stellar effect. In the next section, we will see how to explain some 
unusual afterglows based on the above conclusion. 



3 Some Unusual Afterglows 
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3.1 GRB 980519 



The optical afterglow ~ 8.5 hours after GRB 980519 decayed as oc t" 2 - 05 * - 04 in BVRI [18], 
while the power-law decay index of the X-ray afterglow ax = 2.07 ±0.11 [25], in agreement 
with the optical. The spectrum in optical band alone is well fitted by a power low ^- 12 o±o.25^ 
while the optical and X-ray spectra together can also be fitted by a single power law of the form 
j,-i.05±o.io. i n addition, the radio afterglow of this burst was observed by the VLA at 8.3 GHz, 
and its temporal evolution oc t® 9±0 - 3 between 1998 May 19.8UT and 22.3UT [14]. 

We now analyze the observed afterglow data of GRB 980519 based on our model. We 
assume that for this burst, the forward shock evolved from an ultrarelativistic phase to a non- 
relativistic phase in a dense medium at ~ 8 hr after the burst. So, the detected afterglow, in 
fact, was the radiation from a nonrelativistic shock. This implies 7 ~ 1 at i& ~ 1/3 days. From 
equation (2), therefore, we find 

n 5 ~ 27£ 52 [(l + ^)/2] 3 . (6) 

If p « 2.8, and if the observed optical afterglow was emitted by slow-cooling electrons and the X- 
ray afterglow from fast-cooling electrons, then according to equation (3), the decay index ur = 
(21 — 15p)/10 ~ —2.1 and ax = (4 — 3p)/2 ps —2.2, in excellent agreement with observations. 
Furthermore, the model spectral index at the optical to X-ray band and the decay index at the 
radio band, (5 = —{p — l)/2 « —0.9 and a = 1.1, are quite consistent with the observed ones, 
-1.05 ± 0.10 and 0.9 ± 0.3, respectively. 

We [12] took into account three observed data which correpond to the radio, R-band and 
X-ray frequencies respectively, and inferred intrinsic parameters of the shock and the redshift 
of the burst, e e ~ 0.16, e# ~ 2.8 x 10~ 4 , E 52 ~ 0.27, n 5 ~ 3.4, and z ~ 0.55. After considering 
these reasonable parameters, we numerically studied the trans-relativistic evolution of the shock 
[35] and found that our dense medium model can provide an excellent fit to all the observational 
data of the radio afterglow from GRB 980519 shown in [15]. 



3.2 GRB 000301C 

The optical afterglow data of GRB 000301c were presented in [30]. In addition, the spectral index 
(3 = 1.1 ±0.1. We [13] combined the dense medium model with the pulsar energy injection effect 
to explain the unusual optical afterglow. For stage (ii), if p = 3.4, then a = (12 — 5p)/5 = —1.0 
and (3 = — (p— 1)/2 = —1.2 are consistent with the GRB 000301c R-band afterglow data in initial 
7.5 days after the burst. These data indicate a\ ~ —1-1, which implies a Q b s ~ /?obs at early 
times. If the afterglow were radiated by fast-cooling electrons in the shocked medium, we would 
find a = 2(1 — f3), which is clearly inconsistent with the observational result. Therefore, the 
GRB 000301c R-band afterglow arose from those slow-cooling electrons in the shocked medium. 
For stage (iii), in the case of p = 3.4, the model's time index a = (21 — 15p)/10 = —3.0 is quite 
consistent with the observational data of the GRB 000301c R-band afterglow at late times, 
02 = —3.01 ± 0.53 [30]. We also carried out simulations of the evolution of a shock with energy 
injection from a pulsar and the resulting emission [33]. Our numerical results indeed show one 
sharp break in the late-time afterglow light curve and give a good fit to the R-band afterglow 
data of GRB 000301c. 



5 



4 Conclusions 



We discussed the evolution of an adiabatic shock expanding in a dense medium from an ul- 
trarelativistic phase to a nonrelativistic phase in more details in this paper. In particular, we 
discussed the effects of synchrotron self absorption and energy injection on the afterglow from 
this shock. In a dense medium, the shock becomes nonrelativistic rapidly after a short rela- 
tivistic phase. This transition time varies from several hours to a few days when the medium 
density is from 10 5 to a few xl0 6 cm~ 3 , and the shock energy from 10 51 to 10 54 ergs. The 
afterglow from the shock at the nonrelativistic stage decays more rapidly than at the relativistic 
stage, while the decay index varies from —1.35 to —2.1 if the spectral index of the accelerated 
electron distribution, p = 2.8, and the radiation comes from those slow-cooling electrons. Since 
some models mentioned above predict that a strongly magnetic millisecond pulsar may be born 
during the formation of GRB, we also discuss the effect of such a pulsar on the evolution of 
a nonrelativistic shock through magnetic dipole radiation, in contrast to the case discussed in 
[8,10]. We found that after the energy which the shock obtains from the pulsar is much more 
than the initial energy of the shock, the afterglow decay will flatten significantly and the decay 
index will become —0.4. When the pulsar energy input effect disappears, the index will still be 
—2.1. These features are in excellent agreement with the afterglow of GRB 980519. Further- 
more, our model fits very well all the observational data of GRB 980519 including all the radio 
data. Our model also provides a good fit to the R-band afterglow data of GRB 000301c. 
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